Astronomy & Astrophysics manuscript no. Agudo'efar2006'0716+714 


February 5, 2008 


(DOI: will be inserted by hand later) 





Testing the inverse-Compton catastrophe scenario in the 
intra-day variable blazar S5 0716+71 

II. A search for Intra-day variability at millimetre wavelengths with the IRAM 

30 m telescope 

^ ■ I. Agudo\ T.P. Krichbaum^ H. Ungerechts^, A. Kraus^ A. Witzel^ E. Angelakis^ L. Fuhrmann^'^'"^, U. Bach''"^, S. 
Q ■ Britzen', J. A. Zensus^, S. J. Wagner^, L. Ostorero^'^, E. Ferrero^, J. Gracia^, and M. Grewing^ 

(N 

P5 ' ^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel, 69, D-53121, Bonn, Germany 

^ , ^ Institute de Radio Astronomia Milimetrica, Avenida Divina Pastora, 7, Local 20, E- 18012, Granada, Spain 

Dipartimento di Fisica, Universita di Perugia, via A. Pascoli, 06123 Perugia, Italy 
I INAF, Osservatorio Astronomico di Torino, via Osservatorio 20, 10025 Pino Torinese (TO), Italy 

Landesstemwarte Heidelberg-Konigstuhl, Konigstuhl, D-69117, Heidelberg, Germany 
Tuorla Observatory, University of Turku, Vaisalantie 20, 21500 Piikkio, Finland 
Department of Physics, University of Athens, Panepistimiopolis, 157 84 Zografos, Athens, Greece 

Institut de Radio Astronomie Milimetrique, 300 Rue de la Piscine, Domaine Universitaire de Grenoble, St. Martin d'Heres, 
F-38406, France 



3 
4 
5 
6 

> ■ 

^ 
o 
\o 
o 

Received December 6, 2005; accepted May 30, 2006 

P 

Abstract. We report on a densely time sampled polarimetric flux density monitoring of the BL Lac object 85 0716+71 at 
^ ^. 86 GHz and 229 GHz. The source was observed with the IRAM 30 m telescope at Pico Veleta within a coordinated multi- 

I ' frequency observing campaign, which was centred around a 500 ks INTEGRAL observation during November 10 to 16, 2003. 

O I The aim of this campaign was to search for signatures of inverse-Compton catastrophes through the observation of the broad- 

band variability of the source. At 86 GHz, S5 0716+71 showed no intra-day variability, but showed remarkable inter-day vari- 
^ ' ability with a flux density increase of 34 % during the first four observing days, which can not be explained by source extrinsic 

causes. At this frequency, making use of a new calibration strategy, we reach a relative rms accuracy of the flux density measure- 
ments of 1.2%. Although the flux density variability at 229 GHz was consistent with that at 86 GHz, the larger measurement 
KJf ' errors at 229 GHz do not allow us to detect, with high confidence, inter-day variations at this frequency. At 86 GHz, the linear 

, polarization fraction of S5 0716+71 was unusually large (15.0 ± 1.8) %. Inter-day variability in linear polarization at 86 GHz, 

' with significance level i 95 %; crp/ < P >= 15 % and <r^ = 6°, was observed during the first four observing days. From the 

total flux density variations at the synchrotron turnover frequency (~ 86 GHz) we compute an apparent brightness temperature 
Tg'''' > 1.4 X lO'"* K at a redshift of 0.3, which exceeds by two orders of magnitude the inverse-Compton limit. A relativistic 
correction for T^''' with a Doppler factor S > 7.8 brings the observed brightness temperature down to the inverse Compton 
limit. A more accurate lower limit of <5 > 14.0, consistent with previous estimates from VLBI observations, is obtained from 
the comparison of the 86 GHz synchrotron flux density and the upper limits for the synchrotron self-Compton flux density 
obtained from the INTEGRAL observations. The relativistic beaming of the emission by this high Doppler factor explains the 
non-detection of "catastrophic" inverse-Compton avalanches by INTEGRAL. 

Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: S5 0716+71 
- radio continuum: galaxies - radiation mechanisms: non-thermal - polarization 



1. Introduction 

The effect of intra-day variability (IDV) of radio loud active 
galactic nuclei (AGN) in the radio bands was discovered in 
1986 (Witzel et al. 119861 Heeschen et al. ll987» and since then 
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led to controversial discussions about its physical origin. Most 
of the objects presenting radio IDV appear compact -when are 
observed in the radio and optical bands-, their relativistic jets 
are highly core-dominated on VLBI scales and exhibit high 
brightness temperatures (e.g. Wagner & Witzel 1995i. IDV is 
a common phenomenon in extragalactic flat spectrum radio 
sources and is observed, at centimetre wavelengths, in about 
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10% to 25% of all objects of this class (Quiirenbach et al. ll992l 
Kedziora-Chudczer et al. 120011 Lovell et al. l2003t . The "clas- 
sical IDV" (of type II, as defined in Heeschen et al. 1987 ) is 
characterised by variability amplitudes ^20% and variability 
time-scales between 0.5 days and 2 days. In parallel to the vari- 
ability of the total flux density, similar or even faster variations 
of the linear polarization (Quirrenbach et al. 1 19891 Kraus et al. 
I1999al [r999b 2003 1 and of the circular polarization (Macquart 
et al. l2000t are observed. Intensity and polarization variations 
can be correlated or anti-correlated (e.g. Wagner et al. il996i 
Macquart et aL EUUni Qian et al. EUU^ . 

Observations of IDV sources in the radio and millime- 
tre bands are particularly important, since they reveal the 
highest apparent brightness temperatures {T^ oc y^^; e.g. 
Readhead 1994 1. However, the observed radio emission might 
also be severely affected by interstellar scintillation (ISS), 
which is unavoidable because of the small intrinsic sizes of IDV 
sources (Rickett[l990). It is known that, for high Galactic lati- 
tude blazars, the amplitude of IDV due to the effect of interstel- 
lar scintillation rapidly decreases with frequency for v i 5 GHz 
(e.g. Rickett 1990 Rickett et al. [T9951 Beckert et al. 2002 1. 
In such case and for a source whose size is constant with fre- 
quency, the ISS-induced variability amplitude scales as oc v"^ 
(e.g. Beckert et al. 2003- Hence, any rapid and large amplitude 
source variability seen at millimetre (or shorter) wavelengths 
should not be strongly affected by ISS. Therefore, monitoring 
of IDV sources at mm-bands is particularly important to de- 
termine whether the observed IDV is mainly source extrinsic 
(scintillation), due to source intrinsic processes, or a mixture of 
bofli. 

S5 0716+71 (hereafter 0716+714) is an extremely active 
BL Lac object which varies on time scales from less than one 
hour to months, from radio to X-rays (e.g. Wagner et al. .l996. 
Kraus et al. [2003 Raiteri et al. l2U03t . Alfliough flie optical 
spectrum of 0716+714 appears featureless, even with 4 m class 
telescopes, the absence of any signature of a host galaxy (in 
deep images) sets a lower limit to its redshift of z > 0.3 
(Wagner et al. 1996V Therefore, only lower limits to the bright- 
ness temperature of the source can be derived. The strongest 
constraint on the brightness temperature of this source comes 
from the observed short variability time-scales at radio wave- 
lengths, which sets, via the causality argument, an upper limit 
to the source size of a few tens of micro-arcseconds. During 
the last decade the source has been observed repeatedly in si- 
multaneous centimetre wavelength campaigns (Wagner et al. 
I1990II1996I Otterbein et al. lT^ Kraus et al.|2003 Raiteri et 
al.|2003V during which it typically showed IDV. 0716+714 is 
the only IDV source which has shown simultaneous variability- 
time-scale transitions in the radio and optical bands, which 
have been interpreted as evidence for source intrinsic variabil- 
ity (Quirrenbach et al. 1991, Wagner & Witzel 1995. Qian et 
al. ll996> . The correlation between the radio spectral index and 
the optical flux suggested that the radiation was produced by 
the same particle population (Qian et al. .1996. 1. Additional ar- 
guments for the relevance of the intrinsic origin of the IDV 
in 0716+714 come from the variability amplitude, which in- 
creases from radio to optical, and the recently detected IDV at 



a wavelength of 9 mm (Kraus et al. '20031. Both effects contra- 
dict expectation for the frequency dependence of ISS. 

If IDV is produced by the intrinsic properties of AGN, very 
small angular sizes of the emitting regions (of a few micro- 
arcseconds) are inferred from their short variability time- 
scales. This usually implies apparent brightness temperatures 
several orders of magnitude larger than the inverse-Compton 
(IC) limit ( < lO'^K, Kellermann & Paulinv-Toth [1^691 
Readhead 1994 Kellermann 2003 1. For incoherent synchrotron 
sources of radiation, this limit could be violated only during 
short time ranges, because the high energy photon densities at 
the IC-limit should lead to the rapid cooling of the emitting 
region through the IC scattering of the synchrotron radiation. 
This process is known as the inverse-Compton catastrophe . The 
relativistic beaming of the radiation coming from the emitting 
region (Rees .1966j has been proposed as a likely cause to ex- 
plain the systematic violation of the IC-limit (e.g. Wagner & 
Witzel [T995> . Other scenarios, as coherent synchrotron radia- 
tion, which is an unconventional but not impossible process in 
AGN, have been proposed (e.g. Benford [l992> . However, the 
inverse-Compton catastrophe scenario can not still be ruled out 
(Wagner & Witzel .1995 ) . 

A coordinated broad-band observing campaign -centred 
around a 500ksec INTEGRAL' observation during November 
10 to 16, 2003- was performed to search for signatures of 
inverse-Compton catastrophes on a flaring state of 0716+714 
through the correlation of the broad-band variability of the 
source. The ground-based observations were performed with 
the VLBA, the Effelsberg 100m, the Metsiihovi 13.7 m, the 
IRAM 30 m, the JCMT, the HHT, and the Kitt Peak 12 m tele- 
scopes and the optical-IR telescopes of the WEBT^ collabora- 
tion. The first analysis of the broad-band (from radio to soft 
y-rays) observations (Ostorero et al.' 2006) do not show obvi- 
ous correlation between the intra-day optical variability and the 
inter-day radio variability displayed by the source. Although 
the apparent brightness temperatures of 0716+714 largely ex- 
ceeded the IC-limit at radio wavelengths, no evidence of IC 
avalanches was seen in the INTEGRAL data. 

In this paper we report on the results and implications 
of the first polarimetric millimetre wavelength IDV observa- 
tions of 0716+714, which were performed (within the above- 
mentioned campaign) with the IRAM 30 m radio telescope 
during the period November 10 to 18, 2003. A combined 
discussion of the radio, millimetre and sub-millimetre single 
dish data set will be presented by Fuhrmann et al. l |inprep.^ . 
Forthcoming papers will describe more detailed analysis of the 
optical observations and more sophisticated theoretical mod- 
elling of the broad-band data set. 

2. Observations and data reduction 

2.1. Observations 

The observations presented here were specially planned to 
achieve accurate measurements of the total flux density with a 
typical temporal resolution better than 1 h. To achieve this goal, 

' INTEmational Gamma-Ray Astrophysics Laboratory 
- Whole Earth Blazar Telescope 
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we observed 0716+714 and a number of suitable calibrators 
continuously for ~ 18 h on each day during November 10 to 18, 
2003, using the IRAM 30 m telescope at Pico Veleta (Granada, 
Spain). The observations were interrupted by intervals where 
the elevation of the program source was lower than 20° or when 
the weather conditions were inappropriate for 86 GHz observa- 
tions. Under these conditions, accurate antenna and gain cali- 
brations were not possible. 

The receiver-cabin optic system of the IRAM 30 m tele- 
scope is optimised to observe with up to four different receivers 
at the same time (Wild 1999). We made use of the AlOO, A230, 
BlOO and B230 heterodyne receivers of the telescope, which 
simultaneously observed at both 86 GHz (AlOO and BlOO, at 
A ^ 3.5 mm, where A is the wavelength) and 229 GHz (A230 
and B230, at A ^ 1.3 mm). This standard observing set-up uses 
a grid to divide the incoming signal into two orthogonal lin- 
ear polarizations and then two Martin-Pupplet interferometers 
(one per polarization) to split the 86 GHz and the 229 GHz 
frequency bands. Hence, at each observing band, the A- and 
B-receivers were sensitive to two orthogonal linear polariza- 
tions. This allowed us to obtain linear polarization information 
at 86 GHz. Receiver A230 was not stable enough to provide 
reliable measurements and hence we lost the polarization in- 
formation at 229 GHz. The observing bandwidths at 86 GHz 
and 229 GHz were 0.5 GHz and 1 GHz, and the typical single 
side band system temperatures were ~ lOOK and ~ 400 K, re- 
spectively. 

All the measurements were performed in "beam-switching" 
mode, using a chopper wheel to subtract the emission con- 
tribution from the sky. During the continuous chopping (at 
a frequency of x 7 Hz) the telescope beam was moved to a 
70" off-position. The individual measurements were performed 
by repeated cross-scans over the source position in azimuth 
and elevation. The number of sub-scans in each measurement 
and for each source was chosen according to the flux den- 
sity of the source, with more sub-scans on the fainter objects. 
For technical reasons, the number of sub-scans were multi- 
ples of four and ranged between 4 and 16. 0716+714 was 
observed with a "duty cycle" of about 1 to 2 measurements 
every 60 minutes. Also six bright and point-like extragalac- 
tic radio sources (S5 0212+735, S5 0633+734, S5 0836+710, 
S4 1642+690, S5 1803+784, S5 1928+738, which will be 
called hereafter by their catalogue codes) were observed with a 
duty cycle of at least one measurement every hour Their known 
variability characteristics at cm-wavelengths and location on 
the sky allowed us to assume that they were non-variable on 
time-scales of days and so they were suitable as secondary cal- 
ibrators-'. For the determination of the absolute flux density 
scale, we observed two planets (Mars and Uranus) at least once 
per day and we also included two bright and compact H II re- 
gions (W3 OH, K3-50A) and the planetary nebula NGC 7027 
in the measurement cycle every 4 h to 5 h. 

The frequent observation of the target source and the pri- 
mary and secondary calibrators, in combination with our non- 
standard data reduction procedure, enabled us to calibrate the 

^ In ^3.1l we show that this assumption was true for all the observed 
extragalactic calibrators but 1642+690 



measurements with unprecedented precision. This was possi- 
ble thanks to the accurate determination of the dependence of 
calibration on elevation and time from the measurements of the 
calibrators. This method has been successfully proven on previ- 
ous observations of IDV (e.g. Quirrenbach et al 1992 Kraus et 
al. 120031 . which achieved relative calibration accuracies better 
than 1 % at radio frequencies. The method also allows for the 
empirical determination of the internal calibration errors with- 
out the need of a detailed knowledge of the individual sources 
of error aff'ecting the measurements. 

2.2. Total flux density data reduction 

The initial step of the data calibration used a method equivalent 
to the "chopper-wheel"*" calibration (e.g. Kutner & UIich ll9}?Tl 
to translate the detected power (measured in arbitrary units) 
into calibrated and opacity corrected antenna temperatures (T^ 
in K). The "chopper-wheel" procedure uses measurements of a 
hot and a cold load, both at known temperatures, and the sky 
temperature to compute the calibration factor for each measure- 
ment. At the IRAM 30 m telescope, this calibration factor takes 
into account the total absorption and thermal emission of the 
atmosphere, which are computed by an atmospheric radiative 
transfer model (see Kramer 1997jand references therein for 
details). For our adopted observing set-up at the IRAM 30 m 
telescope, this is the standard initial calibration, which is per- 
formed in "real time" by the software of the telescope. 

The post-observation data reduction was performed follow- 
ing an incremental strategy in different stages, further improv- 
ing the calibration accuracy after each of them. After each cor- 
rection the data were inspected and edited following differ- 
ent criteria, depending on the previous correction. The rms of 
the normalised measurements of the secondary calibrators are 
listed in Table^to show how much each of the corrections im- 
proved the calibration. Here we explain the coiTections applied 
on each of the data reduction stages: 

Stage 1: We first averaged for each scan independently 
the sub-scans in azimuth and in elevation. Then, we fitted a 
Gaussian profile to the average of each of the two scanning di- 
rections. The amplitude of the Gaussian coiTesponds to T* in 
both directions. 

Stage 2: Using the measured FWHM of the telescope beam 
(* 28" at 86 GHz and * 1 1" at 229 GHz), we were able to cor- 
rect the amplitudes measured on the two directions for the ob- 
served (small, typically i 3") telescope pointing errors. These 
corrections were typically smaller than 5 % at 86 GHz and 25 % 
at 229 GHz. After this correction, we averaged the measure- 
ments from the azimuth and elevation directions to produce a 
single measurement per receiver and scan. 

Stage 3: Subsequently, we corrected for elevation depen- 
dent effects on the data, i.e. caused by gravitational deforma- 
tion of the telescope dish. For this purpose, we used the antenna 
gain curves given by Greve et al. ( 1998 1. After this coiTection 

Note that the "chopper-wheel" calibration procedure (which at the 
30 m telescope does not use a chopper wheel) is totally unrelated to the 
"beam-switching" mode (which uses the chopper wheel). This para- 
doxical nomenclature has been maintained for historical reasons. 
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( i 2 % at 86 GHz and i 5 % at 229 GHz), the data did not 
show any significant elevation dependence. 

Stage 4: After that, we removed the most obvious system- 
atic time-dependent variations, i.e. gain decreases (which were 
i 10 % at 86 GHz and < 30 % at 229 GHz) mainly produced by 
abrupt ambient temperature changes. The corresponding gain 
corrections were computed by making use of the densely time- 
sampled measurements of all our secondary calibrators besides 
1642+690, and assuming that they were constant during our 
observations. The latter was confirmed by our variability anal- 
ysis (see § I3.1> . 

Stage 5: At 86 GHz, it was also possible to correct for 
a small relative gain difference of receivers A 100 and BlOO 
(^ 1.035%), detected by comparing good quality measure- 
ments of unpolarized calibrators. After that, the averaging of 
the results from A 100 and BlOO enabled to reduce the un- 
certainties on the determination of the T* measurements at 
86 GHz. An accurate relative calibration of the AlOO and B 100 
gains is also required to measure the linear polarization prop- 
erties of the polarized sources at 86 GHz (see AppendixlAland 
§ I2.3> . Due to the lack of reliable measurements from the A230 
receiver, these data improvements could not be performed for 
the 229 GHz data. 

Stage 6: To further improve our calibration, an additional 
time-dependent correction, similar to the one explained on 
stage 4, was applied to correct for residual systematic time vari- 
ations of the overall gain. The corrected gain variations were 
i 5 % at 86 GHz and < 20 % at 229 GHz on this stage. 

Stage 7: Then we derived the Kelvin-to-Janksy conversion 
factors (Ckj) by comparing the measurements of Mars, Uranus, 
W3 OH, K3-50A and NGC 7027 (which are all of them un- 
polarized) with their assumed absolute surface brightnesses 
(TableOJ. The final values, Ckj(86GHz) = (6.37 + 0.04) Jy/K 
and Ckj(229 GHz) = (9.3 + 0.3) Jy/K are the result of the aver- 
age of the different Kelvin-to-Janksy conversion factors com- 
puted from each of the above calibrators. The errors in Ckj 
take into account the uncertainties of our measurements of the 
primary calibrators, but not the possible inaccuracies on the as- 
sumed values of their surface brightnesses. Ckj(86 GHz) and 
Ckj(229 GHz) were finally applied to all the data to obtain ab- 
solute flux density measurements (S in Jy) for 0716+714 and 
the secondary calibrators. 

Note that the absolute surface brightnesses of Uranus, 
W3 OH, K3-50A and NGC 7027 were calibrated relative to 
those of Mars (Table |2l Griffin & OrtonTWT Kramer!^ i. 
Hence, an additional 5 % absolute calibration error, coming 
from uncertainties in the true martian temperatures (Griffin & 
Orton[l993}, would have to be quadratically added to our ab- 
solute flux density results. However, as this does not affect 
our main relative variability analyses -based on measurements 
scaled to the same Kelvin-to-Janksy factor- the extra 5 % error 
has not been added, unless explicitly indicated for each partic- 
ular absolute flux density measurement. 

Stage 8: Finally, we performed an empirical estimate of 
the internal calibration errors, including those that could not 
be taken into account previously (mainly those affecting the 
gain corrections of data reduction steps 4, 5 and 6). This - 
a posteriori- error estimate was performed by computing the 



Table 1. rms of the normalised measurements of the secondary 
calibrators (but 1642+690) after each of the data reduction 
stages. 



Data reduction stage 


Frequency 


1 


2 
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4 5 


6 


7 


8 
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86 GHz 


12 


11 


5 


4 4 


1.4 


1.2 


1.2 


229 GHz 


120 


100 


30 


19 ... 


16 


16 


16 



rms of the normalized flux density measurements of the as- 
sumed non-variable secondary calibrators (all of them except 
1642+690), which resulted rmsge = 1 .2 % and rms229 - 16% 
for the 86 GHz and 229 GHz data, respectively. These quan- 
tities were finally added in quadrature to the error of each in- 
dividual measurement on stage 7, which were initially com- 
puted from the 1-cr uncertainty estimates of the amplitude of 
the Gaussian profiles and were then propagated through the 
different averages and corrections of the data reduction proce- 
dure. The final relative errors of each individual measurement 
were typically of ~ 2% at 86 GHz and ~ 18% at 229 GHz 
for 0716+714. As indicated above, an additional 5 % factor 
should still be added quadratically to account for the absolute 
flux density errors. 

Note that the small receiver-temperature changes moni- 
tored during the observations did not allow flux density fluc- 
tuations larger than ~ 0.5 % either at 86 GHz or at 229 GHz. 
Hence, the residual rmsge and rms229 statistical uncertainties 
are most likely resulting from the small and rapid atmospheric 
fluctuations between independent measurements, which due to 
the limited time sampling could not be corrected. 

It is worth to stress that the small rmsge - 1 .2 % charac- 
terises the excellent performance of the IRAM 30 m telescope 
and the stability of 86 GHz observations at the telescope site. 
This also demonstrates the ability for future high accuracy IDV 
studies at millimetre wavelengths. 

2.3. Linear polarization data reduction 

To extract the linear polarization information at 86 GHz we 
modelled the response of the two orthogonally polarized lin- 
ear feeds of the receivers AlOO and BlOO to a partially linearly 
polarized source^. In AppendixlAlwe show that, for our observ- 
ing set-up at the IRAM 30 m telescope, these responses fulfil 
the following equations: 

S AWo(i) + S BWo(i) „ 

2 ^S()(i), (1) 

5aioo(0 - Soii) = ^P{i) cos 2 («(/) + ^) (2) 

^ We assumed a negligible degree of circular polarization (pc), 
which is null for the observed Galactic sources and typically < 1 % at 
cm and mm wavelengths for extragalactic AGN (Clemens Thum pri- 
vate communication; Homan & Lister 2006 and references therein). 
In particular, at 86 GHz the measured pc of 0716+714 is < 0.4% 
(Clemens Thum private communication) 
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Table 2. Absolute surface brightnesses (in Jy/beam) of Mars, Uranus, W3 OH, K3-50A and NGC 7027 adopted for the computa- 
tion of the Kelvin-to-Janksy conversion factors of the 30 m telescope at 86 GHz and 229 GHz. 



U.T. days of Nov. 2003 

Source 11 12 13 14 15 16 

86 GHz 

Mars" 149.75 146.88 144.08 141.33 138.65 136.02 

Uranus" 7.06 7.05 7.04 7.03 7.01 7.01 

W3 OH'' 3.94 
K3-50A'' 6.29 
NGC 7027'' 4.71 

229 GHz 

Mars" 691.16 683.33 675.51 666.71 659.92 652.17 

Uranus" 32.88 32.83 32.77 32.72 32.67 32.61 

W3 OH'' 6.48 

K3-50A'' 6.98 

NGC 7027'' 3.68 



"Computed from the expressions given in Kramer H997S, and references therein, which take into account the angular size of the planets. The 
values given for Uranus were computed from an empirical model which used Mars as primary calibrator (Griffin & Orton 1993,i. 
''Obtained from averages of measurements performed during several years - before 1996 - at the IRAM 30 m telescope. These measurements 
were themselves calibrated through observations of planets. The sizes of these calibrators at 1.2 mm are given in Lisenfeld et al. <2000> . 



and 

5bioo(0 -5o(0 = ^P(0cos2(a(0), (3) 

where 5aioo(0 and 5'bi()o(0 are, respectively, the flux densities 
recorded by the receivers A 100 and BlOO at the time of the 
;* measurement, So{i) is the total flux density of the source, 
P{i) is its linearly polarized flux density and □'(/) is the angle 
between the horizontal axis of the telescope's polarization ref- 
erence plane (coincident with the BlOO polarization axis) and 
the direction of polarization of the source at time /. The angle 
a{i) is defined as 

a(i) ^ ^ - X(i) + nd) - , (4) 

where is the astronomical electric vector position angle 
(polarization angle hereafter) of the source (defined in the equa- 
torial coordinate system from north to east), rjii) and e(/) are, 
respectively, the parallactic and elevation angles at the time of 
the fi^ measurement (see Appendix IXI and Thum et al. 2000 J- 
We estimated the magnitudes of P(i) and xiO by performing 
independent fits on the equations (|3 and Q to the data. It was 
necessary to assume that the source did not vary in polarization 
during the time interval of each fit, 5 h t i 12 h. To reduce the 
final uncertainty of the results, the two independent estimates 
of P(i) and xii) (from the independent fits of ^ and (|5}) were 
averaged together. 

Thum et al. f2003^ reported that, for polarimetric measure- 
ments with the A 100 and BlOO receivers of the IRAM 30m 
telescope, the main contribution to the instrumental polariza- 
tion is produced by differences in the relative calibration of 
both receivers. As explained in § 12.21 special care was taken 
to remove these calibration differences /) by correcting for the 
independent small pointing offsets (typically ^ 3") and //) 
by correcting for the gain ratio aff'ecting the A 100 and BlOO 



receivers. Hence, the main component of the instrumental po- 
larization is initially removed from the data through our cali- 
bration procedure. 

To test this and to quantify the amount of spurious polar- 
ization still remaining in the data, we performed fits of the nor- 
malized (to 5'o(0) versions of equations and Q to the data 
of the unpolarized calibrators (Mars, Uranus, W3 OH, K3-50A 
and NGC 7027) over the whole observing time range. The re- 
sult, < p >- (0.2 + 0.9) %, provides a measurement of the 
maximum instrumental linear polarization affecting our polar- 
ization fits, p™" - 1.1%, which we consider negligible for the 
purposes of our data analysis. The fact that p^^" is of similar 
magnitude than rmsgf, indicates that, as for the total flux density 
measurements, our polarization accuracy is most likely limited 
by the statistical uncertainties induced by the small and rapid 
atmospheric fluctuations. 

Such uncertainties are computed from the linear polariza- 
tion fits. In the case of 0716H-714, for time periods f i 12 h, they 
lie in the range 1.2 % to 3.4 % for p and between -0.06 Jy and 
~0.12 Jy for P (see § 13. 2> . The additional error (a 5 % of the 
total flux density) coming from the inaccuracy on the absolute 
flux density scaling, has been added in quadrature wherever it 
was necessary for the final P results; i.e. only for non-relative 
variability results. In such cases, we explicitly indicate that this 
additional error was added. 

The accuracy of our polarization angle estimates is better 
tested on sources with large linear polarization fraction, which 
should show large modulation amplitudes in functions (|2j and 
(|3}, and hence should allow for a more accurate measurement 
of X- Among our observed sources, 0716+714 was the most 
suitable one for such test, since it showed the highest polar- 
ization degree (see §|3} and was better time sampled. The two 
independent fits of functions (|2} and Q to the 0716+714 data, 
over ten diff'erent time ranges of f i 12 h (§ 13 .2> . showed an in- 



6 



I. Agudo et al.: A search for intra-day variability at millimetre wavelengths with the IRAM 30 m telescope 



ternal consistency of the x estimates to a precision better than 
0.5°. This demonstrates the accuracy of equations (|2ji, ^ and 
and the good orthogonality of the signals recorded by re- 
ceivers A 100 and BIOO. However, the statistical uncertainties 
on the determination of fail to be better than A'™," = 1 ■ 1 ° even 
for the polarization fits to the 0716+714 data over the whole 
observing time. Also for 0716+714, these uncertainties range 
between 1 .8° and 9° for the ten time spans with f ^ 12 h (§ 13. 2> . 

Note that the accuracy of absolute polarization angle 
measurements obtained from A 100 and BIOO data is driven 
by the accuracy of the orientation of their polarization axes 
(see Appendix |X| and Thum et al. I2000> . Such orientations 
were calibrated with an accuracy of ~ 1° (Clemens Thum, 
private communication), which is therefore the better absolute 
polarization angle accuracy achievable with the A 100 and 
BIOO receivers. A corresponding error of 1° has been added 
quadratically to the x statistical uncertainties of our measure- 
ments obtained for the whole observing time range, but not 
to those computed for shorter time spans -which are aimed at 
relative variability analyses-. 



3. Results 

In Fig. ^we present plots of the resulting calibrated responses 
of receivers AlOO and BIOO to the 0716+714 emission during 
November 11 to 17. The curves show a remarkable increase 
of ~ 1.5 Jy in 5aioo and 5bioo during November 10 to 14 to- 
gether with a large amplitude sinusoidal modulation, which has 
a characteristic periodicity of 12 h and has a phase diff'erence of 
Trrad between the 5 aioo and 5bioo patterns. The main contribu- 
tion to this modulation is explained by the right-hand terms 
of expressions (|2} and (|3}, which are significative when the 
source is highly linearly polarized. Following (|2j and (|3}, the 
linear polarization flux density of 0716+714 can be estimated 
from the amplitudes of the modulations of the 5aioo and 5bioo 
curves, which lie in the range 0.4 Jy ^ P 0.9 Jy. Hence, the 
curves in Fig.[2represent the first evidence of a high degree of 
linear polarization (10% < p i 20 %) in 0716+714 at 86 GHz. 
As far as we know, such a high linear polarization was never 
observed before in the integrated emission of 0716+714 in any 
radio to mm band. 

In Fig. 12] we show the calibrated response of B230 to the 
0716+714 emission during November 10 to 15. The larger 
measurement uncertainties at 229 GHz, which are typically 
~ 18 %, do not allow us to detect IDV with amplitude lower 
than this level. Therefore, we are able to determine only the 
basic emission properties at this frequency. Since at 229 GHz 
0716+714 does not show the expected 12 h periodic polariza- 
tion modulation, then, within the errors, the response of B230 
is a valid measurement of the source total flux density {8229)- 

3.1. Total flux density 

Accurate measurements of the 86 GHz flux densities (5 86(0) 
were obtained by averaging together 5 aioo(0 and 5bioo(0 (see 
AppendixIXt. The resulting evolution of 586(0 during the time 
range of our observations of 07 16+7 14 is presented in Fig.|3]to- 
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Fig. 1. Calibrated measurements of 5aioo and 5bioo for 
0716+714 vs. time during 2003 November 10 to 16. The si- 
nusoidal pattern is due to source polarization (see the text). 
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Fig. 2. Calibrated measurements of 5b23{) for 0716+714, 
0836+710 and 1928+738 during 2003 November 10 to 15. The 
dashed line on the top inset represents the best linear fit to the 
0716+714 light curve. 



I. Agudo et al.: A search for intra-day variability at millimetre wavelengths with the IRAM 30 m telescope 



7 



5.5 
5.0 



3.5 
3.0 



<SjB> = 4.50 ± 0.52 Jy 



I ^1.20 

\ i 1.10 



1 1 

0836+ 



4 15 16 1 

<Sbb> = 1.49 ± 0.02 Jy 



-=- 1 .00 c^" 
z V 

H 0.90^ 
I <^ 
i 0.80 

^0.70 
7 

1.20 

1.10 

A 

1.00 t| 
V 

0.90 

c 

0.80 
0.70 



to 



1 1 
1928+ 



2.4 
2.2 

2.0 



12 13 14 15 16 17 

<Sb6> = 1.87 ± 0.02 Jy 



L tttll ttt|ttt Httttti 



1.20 




1.10 








1.00 


V 


0.90 


\ 




to 


0.80 




0.70 





11 12 13 14 15 16 17 

J.D. - 2452943.5 (UT days of 2003 November) 

Fig.3. Calibrated 86GHz light curves of 0716+714, 0836+710 
and 1928+738 during 2003 November 10 to 16. The dashed 
line on the top inset represents the best linear fit to the 
0716+714 flux density evolution during 2003 November 10 to 
15. 



gether with the light curves of the non-IDV sources 0836+710 
and 1928+738. The latter demonstrate the accuracy of our cal- 
ibration, which provided flux density measurements charac- 
terised by an nns^d w 1.2%. This represents an unprecedented 
stability for 86 GHz measurements with the IRAM 30 m tele- 
scope. 

The variability of 0716+714 at 86 GHz appears to be dom- 
inated by a monotonous linear increase from < >= (3.71 + 
0.08) Jy during November 10 to 1 1 to < 5 >= (5. 13 +0.86) Jy 
during November 15 (with <> denoting averaging in time). 
During the ~ 5h of data obtained on November 16, the 
source displayed a similar mean flux density of < Sg,(, >= 
(5.12 ± 0.13)Jy. The linear fit of the light curve from days 
10 to 15 (Fig.|3|l gives an increasing rate (Ss6/ < Ss6 >) of 
(8.58 ± 0.14)%/day. The fit performed to the 229 GHz light 
curve gives 5229/ < -5229 (6.2 ± 1.7)%/day, which shows a 
similar trend at both observing frequencies. This suggests that 
the variability at both 86 GHz and 229 GHz originates from the 
same region of the source and the same physical process. 

To characterise the flux density variability of each source 
we used the statistical formulation described in Quirrenbach 



et al. J2000t and Kraus et al. ( I2003t . In particular, we made 
use of the modulation index m, the variability amplitude Y, the 



Table 3. Statistical parameters characterising the variabil- 
ity of 0716+714 and the secondary calibrators during 2003 
November 10 to 16. 



Source 


N 


<S > 


OS 


m 


x] 


Y 






[Jy] 


[Jy] 


[%] 




[%] 




V 


= 86 GHz, mo = 


1.20% 






0212+735 


33 


1.101 


0.016 


1.45 


0.3 




0633+734 


25 


0.942 


0.018 


1.90 


0.4 




0716+714 


74 


4.502 


0.516 


11.46 


55.3 


34.19 


0836+710 


42 


1.486 


0.020 


1.36 


0.1 




1642+690* 


24 


1.044 


0.060 


5.73 


1.6 




1803+784 


32 


1.334 


0.016 


1.19 


0.1 




1928+738 


32 


1.866 


0.024 


1.31 


0.1 






V 


= 229 GHz, mo = 


16.0 % 






0212+735 


15 


0.712 


0.094 


13.14 


0.1 




0633+734 


15 


0.621 


0.085 


13.72 


0.1 




0716+714 


68 


3.506 


0.448 


12.77 


0.5 




0836+710 


31 


0.830 


0.123 


14.77 


0.3 




1642+690* 


8 


0.744 


0.254 


34.17 


0.5 




1803+784 


20 


0.940 


0.235 


25.04 


0.8 




1928+738 


20 


0.904 


0.177 


19.54 


0.4 





'Although this source did not pass our chi-squared test, it displays 
much larger m andx^r than the remaining secondary calibrators, sug- 
gesting a higher variability. 

Note: Y computation was only performed for 0716+714 at 86 GHz, 
the only case for which the variability chi-squared test was passed. 



reduced chi-squared x^- ™d the structure function D(t). The 
corresponding definitions are summarised in AppendixiBl 

The computed statistical parameters for 0716+714 and the 
secondary calibrators are listed in Tablel^for both the 86 GHz 
and the 229 GHz data set. The computation of the variability 
amplitude Y was performed only for those sources that passed 
the chi-squared test (variable with significance level > 99.9 %). 
As can be seen in the Table, only 0716+714 showed such sig- 
nificant variability (at 86 GHz, but not at 229 GHz due to the 
larger uncertainties at that frequency). 

Although 1642+690 formally did not pass this test, it dis- 
played larger values of m and x^- than the remaining calibra- 
tors at both 86 GHz and 229 GHz (see Table O, suggesting a 
weaker variability. In fact, its 86 GHz light curve (frequency 
at which the variability analysis is more accurate) is variable 
at confidence level > 95 %. Based on this evidence, we ex- 
cluded 1642+690 as secondary calibrator from the 86 GHz and 
229 GHz data reduction. 

The structure function D(t) for 0716+714 at 86 GHz is 
shown in Fig. |4] The shape of this function provides the vari- 
ability type classification of the observed sources (Heeschen et 
al. 1987. Kraus et al. 2003) . The lack of pronounced maxima 
and of a saturation point in D(t) on time lags shorter than two 
days allows us to classify 0716+714 as an IDV source of type I 
during the time range of our observations at 86 GHz. The local 
maximum appearing at a time lag of ~ 4 days characterises the 
time-scale of the monotonic increase of the flux density ^ge. In 
agreement with this, and based on the same data set, Fuhrmann 
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November 2003 SF 0716+714 86 GHz 




1 2 3 4 5 



Time Lag (days) 

Fig. 4. D(t) function for 0716+714 at 86GHz. Since the data 
are unevenly sampled, we derived D(t) using an interpolation 
procedure. A minimal interpolation interval of ~ 2.5 h (3 to 
5 times larger than the typical S&e sampling time-range) was 
chosen. To assess the errors caused by the interpolation, D(t) 
was derived twice, first by starting at the beginning of the time 
series and second by proceeding backwards. 



et al. ( |in prep.| l report, from a more sophisticated estimate, a 
time-scale of 3.83;^q |^ days. 

3.2. Polarization 

The 86 GHz polarization estimates were obtained by fitting the 
functions and (|3}, which relate directly the amplitudes and 
phases of such polarization patterns to the linearly polarized 
flux density P and electric vector position angle x of the ob- 
served source, respectively. The fits performed for 0716-1-714 
for i 12h time bins (with a minimum of ten points per fit) are 
shown in Fig.ls] This time binning was chosen as the optimum 
to provide more than one P and x measurement per day and 
the lowest possible uncertainties in their estimate. Polarization 
fits for i 24 h time bins and for the entire time range of the 
observations were also performed for both 0716-1-714 and the 
secondary calibrators (Fig.|6l Table |3 and Table |5jl. 

The fits presented in Fig. |5l confirm the previously men- 
tioned high polarization level of 0716-1-714 at 86 GHz. The av- 
erage polarization properties of the source (during our observa- 
tions) were < p >= (15.0 + 1.8) %, or < P (0.68 + 0.23) Jy, 
and < >= 31.4° + 1.5° (Table|4|l, where the eiTors of < P > 
and < X > take into account those from the absolute total flux 
density and polarization angle uncertainties, respectively (see 
§13. Such a high polarization degree is rather unusual for AGN 
at centimetre wavelengths, for which p > 10% is rarely ob- 
served even during flaring states (e.g. Aller et al. _1999 l and 
typically I % ^ p ^ 6%. Synchrotron self-absorption, usually 
dominating at low centimetre radio frequencies, can reduce p 
(e.g. Pacholczyk ,1970) . However, low values of polarization 
degree are observed at centimetre wavelengths even for opti- 
cally thin sources. Numerical simulations show that, in order 
to explain these low levels of polarization in compact radio 



Table 4. Average 86 GHz polarization properties of 0716-1-714 
and the secondary calibrators during 2003 November 10 to 16. 



Source < P > [Jy] < p > [%] <X >[°] 



0212+735 


0.14+0.06 


12 


±2 


-34+6 


0633+734 


0.08+0.06 


8d 


:3 


77+9 


0716+714 


0.68+0.23 


15.0 


± 1.8 


31.4+1.5 


0836+710 


0.08+0.08 


5.3 d 


: 1.6 


-76+9 


1642+690 


04+" °^ 


4d 


:4 


66+18 


1803+784 


0.11+0.07 


8d 


:2 


-76+9 


1928+738 


02+° 
"•"^-0.02 


1.1 


+1.6 
-1.1 





Note: These results have been obtained, for each source, from the po- 
larization fits over the entire time range of the observations. The errors 
were originally computed as 1-tr uncertainties of the polarization fits. 
The errors of the < P > and < x > measurements include those from 
the absolute total flux density and polarization angle uncertainties, re- 
spectively (see §|5J. 

sources, a large fraction of the magnetic field (typically 50 %- 
70 %) should be disordered in direction (e.g. Gomez, Alberdi 
& Marcaide 1993). The 15% polarization degree detected in 
0716+714 might indicate a high level of magnetic field align- 
ment in the region from which the 86 GHz emission originated. 

In addition, one of the six AGN calibrators observed dur- 
ing our observations, 0212+735, also displays a polarization 
degree p > 10% and 0633+734 and 1803+784 show p ^ B, %. 
These values are still high when compared with the typically 
observed polarization degrees at centimetre wavelengths. The 
fact that, at 86 GHz, more than 50 % of the observed AGN dis- 
play larger polarization degrees than those typically measured 
at centimetre wavelengths might suggest that the polarization 
degree in AGN increases with frequency. This would indicate 
higher ordering of the magnetic fields in the deeper regions of 
their jets observed in the millimetre bands respect to centime- 
tre wavelengths. However, it is clear that our sampled set of 
AGN is not statistically significant enough to strongly support 
this hypothesis. The millimetre polarization survey performed 
by Thum et al. (in prep.) over more than 70 AGN during the 
summer of 2006 will certainly improve the statistics. 

Inspection of Fig. |5] indicates time-dependent changes in 
both the amplitude and the phase of the fitted functions, sug- 
gesting that the linear polarization of 0716+714 varied along 
our observations. The P estimates, for both 0716+714 and the 
secondary calibrators, are presented as a function of time in 
Fig. while the x evolution for 0716+714 only is shown in 
Fig- HI Due to the much lower polarization degree of the cali- 
brators, their polarization properties were optimally fitted with 
time bins of < 24 h. This provided one polarization measure- 
ment per day to test the stability of the secondary calibrators 
and it also decreased the fitting eiTors compared to shorter time- 
binnings. To allow for a better comparison with the calibrators, 
the 0716+714 data were also fitted using a < 24h time binning 
(Fig.|7lFig.|H]andTable|5}. 

While the secondary calibrators remained stable in polar- 
ization within the eiTors, 07 16+7 14 displayed polarization vari- 
ability during November 1 1 and 14. To quantify this variability, 
we performed a statistical analysis similar to the one presented 
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Fig. 5. Plots of the fits of functions (|2j and (|3} to the 86 GHz data of 0716+714. Each plot represents a time bin i 12h of 
observing time (labelled within each of them) with a minimum of ten data points. The time spans, ordered from top to bottom 
and from left to right, cover the whole observing time. The whole data set has been plotted together in the right-bottom plot. The 
sinusoidal lines represent the independent non-linear least squares fits of functions and Q. The resulting averaged P and^ 
variables are given within each corresponding plot. Errors were originallycomputed as l-cr uncertainties from the fits. For the 
right-bottom plot, the eiTors from the absolute total flux density and polarization angle uncertainties (see §|3 have been taken int 
account (added in quadrature). 



in ^3.11 Due to the much weaker polarization of the secondary 
calibrators compared to 07 16+7 14, the errors of their P and es- 
pecially their estimates are much larger than for 0716+714. 
Hence, in principle, it would not be reasonable to calculate nip, 
m^, Yp and since they do not take into account the uncertain- 
ties on the measurements and they would not provide reliable 
variabiUty results. However, the mean, the standard deviation 
and the reduced chi-squares of both the P and y distributions 
can still provide useful statistical information. These magni- 



tudes are summarised in Table|5] None of the sources analysed 
passed the chi-squared test neither for P nor for x- However, 
0716+714 displayed relatively large values of x^p = 2.4 and 
X^j^ = 2.2, with corresponding variability significance levels 
> 95 %. Hence, for the particular case of 0716+714 and given 
the relatively large significance of its polarization variability, it 
makes sense to compute both mp - 15 % and = cr^ - 6°, 
which are attributed to the source behaviour during the first four 
observing days (Figs.0and|8j. 
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Fig. 6. 86 GHz polarization patterns of the secondary calibrators. For the sake of clarity, a shift of ;rrad was added to function 
Q to match function (|2}, i.e. 5'gjQQ(0 - 5'()(0 - I /2P{i) cos2{a(i) + njl) (see § I2.3> . For each source, the data from the whole 
observing time-range is presented. The sinusoidal lines represent the non-linear least squares fits to the data. The resulting 
averaged P and x variables are given within each corresponding plot. The errors presented on each plot were computed by 
adding in quadrature the 1-cr uncertainties from the fits and the errors from the absolute total flux density and polarization angle 
uncertainties (see §0. 
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Fig. 7. 86 GHz P curves for 0716-H714, 0836-H710 and 
1803H-784. Time bins of < 24h were used for all the sources. 
For 0716H-714 the results for time bins i 12h are also pre- 
sented. Horizontal bars symbolise the time range over which 
each estimate of P was computed. 



The amount of P and x variability of 0716-1-714 between 
each pair of adjacent measurements of these variables was 
characterised through their fractional variability amplitudes: 



V/'(Af,-i,,) 



2(rP(M) - rP(M_l)) 

P(Af,) + P(Af,-i) 



(5) 



X vs. time for 0716+714 



O O <l2h time binning 
• • <24li time binning 



1T 




11 12 13 14 15 16 

J.D.-2452943.5 (UT days of 2003 November) 

Fig. 8. 86 GHz polarization angle vs. time for 0716-1-714. The 
results obtained for time bins of both < 24h and i 12 h are 
presented. Horizontal bars symbolise the time range over which 
each X estimate was computed. For the sake of clarity, the x 
time series of the calibrators (which are all constant within their 
corresponding errors) have not been plotted due to their large 
uncertainties, typically of 15°-20°. 



and 



Vy(Af,--i,/) 



2{x{^td-x(M_,)) 



(6) 



XiMi) + xi^ti^i) 
where Af, (with / - 1,2, ...,10 for the 0716+714 polariza- 
tion measurements performed with a time binning i 12h) 
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0716+714 polarization fractional variability amplitude vs. time 
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Fig. 9. Amplitude of fractional polarization variability, as a 
function of time. The vertical bars represent the Vp and 
propagated errors whereas the horizontal bars symbolise the 
maximum possible time range over which the variability took 
place. See the text for definitions. 

denotes for the interval over which each one of the polariza- 
tion fits was performed and Af,_i ,■ symbolises the time inter- 
val between each pair of adjacent polarization measurements. 
We considered that an intra-day polarization variation event 
is significant when either Vp or V^, is larger than its corre- 
sponding uncertainties computed by propagating the errors of 
P and X- There is only one time interval (Afej) for which 
this condition is fulfilled. Afg 7 is a time interval of ~ 24 h on 
November 13-14, for which both Vp(At(,j) = (22 + 1 1) % and 
Vj^(At(,j) - (26 ± 12) % show a 2cr variation in P and x, re- 
spectively (see Fig. O. Hence, this simple analysis provided 
indications of IDV on the polarized emission of 0716-1-714 dur- 
ing our observations. Nevertheless, this statement should not be 
taken as a strong evidence of polarization IDV, as no 3cr IDV 
event has been detected. A future dedicated millimetre mon- 
itoring campaign, optimised for high time sampling and high 
polarization sensitivity, will be necessary to confirm and clarify 
whether 0716+714 (or any other AGN) is variable in polariza- 
tion on time-scales of i 24 h. 

4. Discussion 

4.1. Spectrum 

The millimetre data presented here contributed to the assem- 
bling of the broad-band spectrum (from radio to y-rays) of 
0716+714 during November 10 to 16, 2003 presented by 
Ostorero et al. ( 2006 1. The spectrum was inverted between 
the cm- and mm-bands. Its turnover frequency (v,„, at which 
the synchrotron spectrum peaks) was located at v,„ « 86 GHz. 
Between 86 GHz and 667 GHz the spectrum can be fitted by a 
single power law with spectral index a ^ -0.3 (where Sy o: v" 
and V denotes the observing frequencies within this spectral 
range). 

At 86 GHz, the total flux density of the source ( ^3.U ranged 
from Sm = 3.55 + 0.19 Jy at the 10.88 U.T. days of November 



Table 5. Statistical parameters characterising the 86 GHz polar- 
ization variability of 0716+714 and the secondary calibrators, 
but 1928+738, during 2003 November 10-16. 



Source 


M 

JN 


< P > 

[Jy] 


[Jy] 


2 


<x > 

n 


n 


2 


0212+735 


5 


0.122 


0.024 


0.3 


-33.1 


6.0 


0.2 


0633+734 


5 


0.092 


0.021 


0.1 


69.3 


8.5 


0.3 


0716+714° 


10 


0.648 


0.103 


1.8 


30.9 


4.3 


1.7 


0716+714* 


6 


0.639 


0.095 


24 


31.0 


3.3 


2.2 


0836+710 


6 


0.085 


0.020 


0.2 


-76.3 


8.5 


0.1 


1642+690 


5 


0.064 


0.023 


0.1 


54.6 


17.6 


0.3 


1803+784 


5 


0.111 


0.022 


0.1 


76.4 


4.5 


1.2 



"0716+714 results for i 12h. 

''0716+714 results for < 24h are also presented to allow direct com- 
parison with the secondary calibrators. 

Note 1 : < P > and < ^ > in this Table, which are based on P and x 
fittings over < 24 h time bins, are less accurate than those from Table 
El which where fitted over the whole time span of the observations. 
Note 2: The weak linearly polarized flux density of 1928+738 (see 
Table|3 does not allow us to analyse its variability. 

2003 to S,„ = 5.23+0.07 Jy at the 15.00 U.T days of November 
2003. The larger total flux density of the source at 86 GHz 
compared to that at 229 GHz, (Figs. |2] and |3} shows that the 
source synchrotron spectrum was optically thin between these 
two frequencies during the whole time range of our observa- 
tions. To attempt to characterise the source spectral evolu- 
tion, we have computed the spectral index for each pair of 
5 86(0 and ^229(0. where / is the index defining each one of 
the simultaneous 86 GHz and 229 GHz measurements. The re- 
sulting Q'86-229 evolution (shown in Fig.llOt is consistent with 
a decreasing pattern (spectral softening) with constant rate 
Q'86-229 = (-0.03+0.02) d^' and average Q'86-229 = -0.23+0.10 
from the 10.88 to the 15.00 U.T days of November 2003. 
However, the computed Q'86-229 is not statistically significant 
and the possible spectral change of 0716+714 needs to be 
tested through more accurate spectral index data. 

4.2. Intrinsic variability 

Given the known frequency dependence for weak ISS, it is 
rather unlikely that the F = 34 % amplitude inter-day vari- 
ability reported at 86 GHz is caused by this extrinsic effect. 
As outlined in ^ such large variability in the mm-bands can 
not be reproduced by the standard ISS models (e.g. Beckert et 
al.l2()02J. Fuhrmann et al. (in prep. 1 perform a more detailed 
variability study combining all the available simultaneous data 
from 5 GHz to 667 GHz, and show that the inter-day variations 
shown in this paper contradict the predictions of such models. 
We should hence attribute the large amplitude 86 GHz variabil- 
ity of 0716+714 to intrinsic causes. 

Hence, the high amplitude inter-day variability reported 
by Ostorero et al. ( 2006) at 32 GHz and 37 GHz, which well 
matches with the millimetre behaviour shown in this paper, 
was most likely produced in the source itself and not by ISS. 
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07' 5 + 71 4: 86-??9 GHz spectral index evolLition 




11 12 13 14 '5 16 17 

J.D.-2452943.5 (JT days of 2003 November) 

Fig. 10. Evolution of the spectral index age-iag in 0716+714. 
The absolute total flux density calibration errors of 5 86 and 
5 229 (see §13 were taken into account to compute those of 
0^86-229 ■ The dashed line represents the best linear fit to the 
data, which is consistent with a pattern with a rate Q'86-229 = 
(-0.03 + 0.02) d-i and average Q'86-229 = -0.23 + 0.10. 



Ostorero et al. i 2006> also report on the rapid IDV in the opti- 
cal range observed during the time range of our IRAM 30 m 
observations. Although 0716+714 did not show evidence of 
simultaneous radio or mm and optical IDV-time-scale transi- 
tions, as during previous monitoring campaigns (Quirrenbach 
et al. I1991il . this can not be taken as an argument in favour of 
the extrinsic origin of the variability. It indicates only that, dur- 
ing our observations, the radio-mm emission was radiated from 
a diff'erent region, or through a different mechanism, than that 
producing the optical emission. 

4.3. Brightness temperature 

Under the assumption that the main component of the mm- 
variability is intrinsic to the source and via the causality argu- 
ment, it is possible to estimate some of the physical conditions 
governing the source. This can be performed by modelling it as 
a relativistically moving homogeneous sphere (e.g. Marscher et 
al.Q979 and Ghisellini et al. ll993> . Under these assumptions, 
the intrinsic brightness temperature of the emitting region at 
the turnover frequency v,„ (at which the spectral index is zero, 
by definition) can be computed from its variability time scale 
as: 



Tb = 1.77 X 10' 



•612 



l+z 



K, 



(7) 



where S,„ is the flux density of the source at the turnover fre- 
quency in Jy and v,„ is in GHz. z is the redshift of the source, 6 
is an equivalent Doppler factor and 



(8) 



6, i 3.56 X lO''^t,dL-\l+z)6mas 



is an estimate of its angular diameter (Marscher et al. 1 19791 1 
from light-travel time arguments in the absence of coherence, 
fv =1 dln5v/df r'= (< Sy > /ASy)At is the variability time 



scale in days and di^ is the source luminosity distance in Gpc. 
In this paper, we adopt z - 0.3 (the lower limit for the redshift 
of 0716+714, see e.g. |T]and Wagner et aL fTOgt . a Hubble 
constant Hq - 72kms ' Mpc ' and a Friedmann- Robertson- 
Walker cosmology with Q,„ = 0.3 and Qa = 0.7. Within these 
assumptions, d]^ - 1.51 Gpc. 

Our 86 GHz observations provide us the values of 5 m - 
5.23Jy, < S„ >= 4.50Jy, A5„, = 1.68Jy, Af = 4.12days 
and ty = 11.1 days. The only remaining variable in expres- 
sion which can not be directly observed is S. By assum- 
ing a non-moving source, expression Q gives us the apparent 
brightness temperature T'^^ = Tb{6 = 1) > 1.4 x 10"*K at 
z - 0.3, which exceeds by more than two orders of magni- 
tude the IC-limit P™^ ^ 10'^ K (Kellermann & Pauliny-Toth 
[l969l. Hence, 7"^"^-, sets a lower limit to the Doppler factor of 
the emitting region, which should then be > 5.2. Note that 
a more accurate limit for T'g'j'p (of ~ 3 x 10" K) was computed 
by Readhead 1994 (see also Kellermann I2003> . To avoid the 
violation of the latter, 6tj, > 7.8 is required. 



4.4. Angular size 

The latter constraint {6tb > 7.8), together with expression 
allows us to estimate an upper limit for the size of the emitting 
region responsible for the variability during our observations, 
which should be Otb < 0.03 mas. Obviously, flj-g is also an up- 
per limit for the maximum size of such region obtained from 
St, > 5.2. 

It should be stressed that 0Tg does not account for the size 
of the whole source, but only for the region responsible for 
the rapid variability, which in the case of 0716+714 has been 
associated with the VLBI core component of its relativistic jet 
(Bach et al. l 2006> . 

An independent constraint on the size of the variability 
region can be obtained from the 86 GHz VLBI observations of 
0716+714 performed by Bach et al. ( 2006) in April 2003. A 
Gaussian model fit of their resulting image gives a FWHM size 
of the 86 GHz VLBI core of 0.014 mas, which is much lower 
than the minimum beam size of VLBI observations at such 
frequency (~ 0.04 mas; e.g. Agudo et al. .2005 J . This gives us 
GvLBi < 0.04 mas, which is in agreement with our estimate of 

We will hence assume hereafter that the true angular size of 
the rapid variability region is 9y = 9t^ < 0.03 mas. 

4.5. Inverse Compton Doppler factor 

It is also possible to compare the expected and the observed 
first order self-Compton y-ray flux density to constrain the 
IC Doppler factor (die) of the source (Marscher .1983. and 
Ghisellini et al. ll993> . die is defined as: 



3IC 



ln(V2/Vm)Vy" 



1/(4-20-) 



(l+Z): 



(9) 



where a is the spectral index of the optically thin synchrotron 
spectrum {a - -0.3 has been adopted here), f{a) - 0.14 - 
0.08a (Ghisellini et al. .1993i and references therein), V2 is the 
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Table 6. Physical parameter estimates. See the text for defini- 
tions. 



Parameter 



Value 



Assumptions 



Forz = 0.3,4 = 1510Mpc, S„, = 5.23 Jy, U. = 11.1 days 



^ B 


> 1.4 X 10"*K 


S = 1 




6Te 


> 5.2 


^bTc = 10'^ K 




Stb 


> 7.8 


r^™c = 3xl0"K 




Stb 


< 0.03 mas 


' (5>7.8 






< 0.04 mas 


86 GHz VLBl beam-size 




e. 


< 0.03 mas 






<5lC,8keV 


> 15.8 


Sgkev < 1-37 X 10^' Jy, a = - 


-0.3 


<5lC,23keV 


> 14.6 


523kev < 2.55 X lO"' Jy, a = ■ 


-0.3 


<5lC,63kcV 


> 15.1 


563kev < 1-23 X lO'' Jy, a = ■ 


-0.3 


<5lC,141kcV 


> 14.0 


5i4ikev < 2.42 X lO"' Jy, a = 


-0.3 



upper frequency cut-ofF of the synchrotron spectrum in GHz, 
Vy is the y-ray energy in keV and Sy is the measured y-ray 
flux at Vy in Jy. V2, in GHz, is the most poorly known vari- 
able among those in equation (|9jl. However, this problem is 
overcome by the fact that 6jc is a weak function of V2, which 
is contained within a logarithmic factor Here, we have as- 
sumed V2 = 10*'GHz (Biermann & Strittmatter 1987). Implicit 
to equation (|5Jl is the assumption that the electrons radiating 
the synchrotron emission have a power-law energy distribution 
(N{E) - NqE^"'^, with hv„, ^Et hv2) and that the magnetic 
field is tangled in the rest frame of the emitting plasma. 

By taking into account the dependence of 0i, on 6 in 6ic 
can be written as: 



(5ic I 



'\n{v2/v„,)vy"[f(a)SM+z)f 



(4-2<r)^l/(10-6a) 



Sy@, 



(6-4ff)y(5-3a) 



(10) 



where we have defined 0,, = 3.56 x lO-\di^-Ul + z) mas. 

The results of the simultaneous INTEGRAL soft y-ray ob- 
servations of 0716+714 reported by Ostorero et al. ( 2006 1 

{F3-35keV < 6.12 X lO^'^ergCm^-S^', Fi5_40keV < 1-41 X 



10"" ergcm"^ s~', F4o-iookeV < 1.72 x 10"" ergcm"^ s"' and 
fioo-200keV < 5.73x 10"" erg cm"^ s"' ; see also Table|6lfor the 
corresponding flux densities in Jy), together with expression 
dlOL allow us to set new constraints on the Doppler factor of 
0716+714 during our observations. Assuming f,, - 11.1 days 
as in §|0] we derive lower limits of 5ic,8kev > 15.8, 5ic,23kev > 
14.6, 5ic,63kev > 15.1 and 5ic,i4ikev > 14.0 from the upper lim- 
its of the flux densities of the source at 8 keV, 23 keV, 63 keV 
and 141 keV, respectively. 

It is worth to note that, for a - -0.3, expression (I10> is 
a weak function of the redshift and the luminosity distance of 
the source (Sic oc (1 + zy^'^djj''^). Hence, the constraint on 
the Doppler factor of the source provided by ( I10> . Sic > 14.0, 
is more robust than those derived from expression (0 and the 
limits for^B (§0- 

5. Conclusions 

We have presented the results from millimetre observations of 
0716+714, which were performed on 2003 November 10 to 18 



with the IRAM 30 m telescope. Our observation strategy, based 
on the rapid time sampling of both the target source and the 
calibrators, enabled us to reach a relative calibration accuracy 
of 1 .2 % at 86 GHz, which demonstrates the good performance 
of this telescope and its ability for future accurate IDV studies 
in the millimetre range. 

During our first four observing days, the source displayed 
large amplitude (Y - 34%) and monotonous inter-day vari- 
ability at 86 GHz. As such large amplitude variability is not 
expected to be produced by standard ISS at mm wavelengths 
(Rickett et al. II995> . this variation should be considered as 
intrinsic to the source and not due to the influence of the in- 
terstellar medium. The similar 32 GHz and 37 GHz behaviour 
reported by Ostorero et al. ( 2006 1 during the same observing 
time range could hence be explained by intrinsic causes also. 

Ostorero et al. i 2006j also report clear evidence of IDV in 
the optical range, which is not matched either in the mm or in 
the cm bands (see also Fuhrmann et al. in prep.^ . This indicates 
that the radio-mm emitting region was located in a different re- 
gion than the optical one or that their radiation behaviours were 
driven by different physical processes during our observations. 

We have reported an unusually large linear polarization 
degree < p >= (15.0 + 1.8)% of 0716+714 at 86 GHz, which 
suggests a large level of magnetic field alignment. At such 
frequency, linear polarization inter-day variability, with signif- 
icance level Z 95 %; crp/ < P >- 15 % and cr^ - 6°, was 
observed during the first four observing days. We have also 
shown 2cr evidence of simultaneous polarization flux density 
and polarization angle IDV within a time i 24 h. If such rapid 
polarization variations that are uncoiTelated with the total flux 
density variability are confirmed by future observations, then 
equally rapid changes of the magnetic field configuration of the 
source or changes of opacity around v„, would be required to 
explain the phenomenon. In both cases, inhomogeneous mod- 
els would probably have to be invoked. 

The synchrotron spectrum of 0716+714 peaked at v,„ x 
86 GHz during 2003 November 10 to 15 (Ostorero et 
al. 1 2006> . and showed an optically thin spectral index 
0-86-229 = -0.23 + 0.10 between 86 GHz and 229 GHz. The ap- 
parent brightness temperature derived from our 86 GHz light 
curve, TgPP > 1.4 x lO'^^K for a redshift z = 0.3, ex- 
ceeds, at least by two orders of magnitude, the IC-limit of 
3 x 10" K (Kellemann & Paulinv-Toth .1969. Readhead J^^?! 
Kellermann l2003t . This mismatch can be explained by the rel- 
ativistic motion or expansion of the source with a minimum 
Doppler factor 6 > 7.8. The upper limits from the soft y-ray si- 
multaneous INTEGRAL observations in the 3 ke V to 200 keV 
energy range enabled us to compute an independent and more 
robust limit for the source Doppler factor, die i 14. This limit 
is consistent with previous estimates of the Doppler factor 
of 0716+714 measured from the kinematics of VLBI-scale jet 
features (Bach et al. 2005 1, which ranged from 5vlbi = 20 to 
i5vLBi = 30. Such high Doppler factors have also been mea- 
sured with 43 GHz-VLBI by Jorstad et al. J2005> . who detected 
^VLBi > 20 in 9 of the 13 monitored blazars. 

As no soft-y-ray IC-avalanches were detected by 
INTEGRAL during our observations and the reported 
large amplitude 86 GHz variability can not be ascribed to the 
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ISS, the relativistic beaming of the radiation coming from 
the emitting region in 0716+714 offers a robust explanation 
of the apparent violation of the IC-limit of the brightness 
temperature in the mm range. Note, however, that (total or 
partial) coherent synchrotron-emission scenarios can not be 
ruled out. 

Finally, we should stress that we have proven that 
07 16+7 14, in particular, and blazars, in general, can display ap- 
parent brightness temperatures two orders of magnitude larger 
than the theoretical limits in the millimetre range and that 
the influence of the interstellar medium is not always neces- 
sary to explain the mismatch between observations and theory. 
Following the above arguments, we are rather confident that, 
apart from inaccuracies in our assumptions, the estimates and 
limits for the physical parameters of 0716+714 reported in the 
previous sections correspond to those governing the observed 
source behaviour 
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Appendix A: Response of the IRAM 30 m 
heterodyne receivers to a partially linearly 
polarized wave 

The power received by an antenna when it observes a source of 
arbitrary polarization can be characterised as 

W ^ -SoAeil - d) + SoAed cos- ^, (A.l) 

where So is the total flux density of the source, Ag is the effec- 
Xive aperture of the antenna, d the source degree of polarization 
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and y the angle between the wave and the antenna polarization 
states on the Poincare sphere (see e.g. Kraus 1986 l. The first 
and second terms on the right-hand side of Eq. jA.U represent 
the unpolarized and polarized responses of the antenna, respec- 
tively. 

If the polarized component of the incident radiation is only 
linearly polarized, then d is the degree of linear polarization 
(p). In addition, for the case of the IRAM 30 m millimetre ra- 
dio telescope, when heterodyne receivers are used, y = a + /3. 
Here, a is the angle between the polarization direction of the 
incident signal and the horizontal direction in the telescope re- 
ceiver cabin -the x axis of the reference polarization plane, see 
Fig. lA.U . y6 is the angle between the x axis and the polarization 
orientation of the receiver For the particular case of receivers 
A 100 and BlOO, /? is k/2 and 0, respectively. 

Taking into account that W = kT^ and T^S obs = Ae/{2k), 
where k is the Boltzmann constant, is the calibrated antenna 
temperature outlined on ^2.2l and S obs is an equivalent total flux 
density observed for the receiver, Eq. dA.H may be expressed 
forAlOO and B 100 as: 

5Aioo(f) = 5o(f)(l - pit)) + 2So(t)p(t) sm\a(t)) (A.2) 
and 

SBiooit) = Soim - pit)) + 2Soit)pit) cos\ait)) . (A.3) 

where we have introduced the possible time dependence of 
each variable and we have defined 5ai(K)(0 = S obs{o'{t),/3 = 
7T/2, t) and 5Bi{)()(f) = S„t.Xa{t),/3 = 0, t). From dX^ and 
and defining the linearly polarized flux density of the observed 
source as P = S op, then it is straightforward to show that: 



zenith 



I AlOO' 



(f) + 5B100(f) 



So(t), 
L 

SAioo(t) - Soit) = P(f)cos2(a(f) + ^) 
and 

5Bioo(f) - 5o(f) = Pit) cos 2(a(t)) , 



(A.4) 
(A.5) 

(A.6) 




E-vector 



- X 



R.A. 



Fig.A.l. Angles in the reference polarization plane (x, y) as 
seen from the AlOO and BlOO, IRAM 30m heterodyne re- 
ceivers. a(t), which is defined in (x, y) and can be measured 
by fitting equations iA.5\ and iA.6\ . relates to the astronomical 
electric vector position angle xit) (defined from north to east in 
the equatorial coordinate system (R.A.,Dec.)), to the parallactic 
angle 77(f) and to the elevation angle e(f). 

Appendix B: Statistical formulation for variability 

The statistical variables used in this paper are defined in detail 

Here we 



(B.l) 



by Quirrenbach et al. J2000> and Kraus et al. J2003> 

present a summary of these definitions. 
The modulation index. 



m[%] = 100 



<S >' 



(where os is the standard deviation of the flux density measure- 
ments and < 5 > is its average in time) provides the strength of 
the observed variability amplitude of each source without tak- 
ing into account the uncertainty of the measurements. 
The variability amplitude. 



Thum et al. I2OOOI showed that, for the IRAM 30 m tele- 
scope, a(t) can be expressed as the following function of astro- Y[%] - 3 Irn^ - 1 
nomical angles: 



(B.2) 



a(t) = 2 - X(t) + riit) - e(t) , 



(A.7) 



(see also Fig. lA.U where xi^) is the electric vector position 
angle of the source (defined from north to east in the equato- 
rial coordinate system (R.A. , Dec.)) and 77(f) and e(f) are, re- 
spectively, the parallactic angle and elevation angles (see Fig. 
lA.ll and. Thum et al. 2000 1. Note that the location of the re- 
ceiver cabin in the Nasmyth focus of the IRAM 30 m telescope 
causes a rotation of the equatorial reference system by an an- 
gle xoit) - e(0 - '7(0 with regard to the polarization reference 
system (x, y). 



(where mo is the modulation index of the non-variable and 
nearby calibrators) gives a measurement of the variability am- 
plitude with regard to that of the non-variable calibrators, for 
which Y is set to zero. 

The reduced chi-squared. 



1 ^ I 



Si -<S > 
ASi 



(B.3) 



(where is the number of measurements, the Si are the in- 
dividual flux densities and A, their errors) tests the hypothesis 
that the measured light curves could be modelled by a constant 
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function. We consider that a source is variable (with a signifi- 
cance level > 99.9 %) only if its light curve shows a probability 
to be constant < 0.1%. 

Finally, the structure function is defined as 

DiT) = <iSit)-Sit + T)f>, (B.4) 

(with <> denoting averaging in time). A characteristic time- 
scale in the light curves (i.e., the time between a maximum and 
a minimum or vice-versa) causes a maximum in D{t), while a 
periodic pattern produces a minimum. 



